Dynamics of vortex matter in rotating two-species Bose-Einstein condensates 
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In a rotating two- component Bose-Einstein condensate (BEC), the traditional triangular vortex 
lattice can be replaced by a rectangular vortex lattice or even a structure characterized in terms 
of vortex sheets, depending on the interspecies interactions. We study the dynamics of this system 
by analyzing the Bogoliubov excitation spectrum. Excitations familiar to BEC vortex systems are 
found such as Tkachenko modes, hydrodynamic modes and surface waves, however, the complex 
two-component morphology also gives rise to new phenomena including shear flow between vortex 
sheets. 
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Since the realization of vortices in trapped atomic 
Bose-Einstein condensate (BEC) 0,0,0, there have been 
extensive efforts to understand the dynamics of quantized 
vortices in a dilute bosonic system from both the theo- 
retical and experimental view 0, . With the exciting 
developments in multi-species BEC, it is natural to ask 
about the nature of vortices in these systems. In the 
absence of rotation, it is known that the ground state 
of a two-species condensate goes through phase transi- 
tion from a phase-overlapped (PO) to a phase-separated 
(PS) condensate depending on the interspecies interac- 
tion strengths 0, Q ;the equivalent transition for a rotat- 
ing system is that from a rectangular lattice to vortex 
sheets The rectangular vortex lattice has been ob- 
served experimentally at JILA with rubidium hyperfine 
states as the two species 9]. The vortex-sheet configu- 
ration was also observed in a different context, ^He-A, 
in which two different spin states of the Cooper pairs 
constitute the two species 0. In this paper we discuss 
the excitations of the rotating two-species condensate in 
both the PO and PS regimes. 

With two mean field order parameters in a 
single system, the inter-species interaction energy, 
gi2 J dr ij)2^\'ipi'ip2 is introduced, where gi2 is the in- 
teraction strength between species 1 and 2. We ignore 
any coherent coupling between species in this paper. 
With this additional term in the free energy, the Gross- 
Pitaevskii equation (CPE) in a rotating frame [TJ be- 
comes (i = 1, 2) 

{n + g^\^^^f + gi2\^^3-if - ^C,) = (1) 

which can be solved to find the ground state of a con- 
densate. Here, H = -{fi^ /2M)V'^ + Vtr(r) and = 
—ifid/dcj). Vt is the rotational frequency of the reference 
frame and gi = ^iraih^ / Mi is the interparticle interaction 
strength of species i with the 5-wave scattering length a^. 
With this choice of parameters, in the absence of rota- 
tion, the transition point from a PO to PS condensate is 
found at gi2 = g 0, 13- 

The coupled Bogoliubov-de Gennes (BdG) equations 
[m that describe linearized dynamics around the ground 



states ipi and iIj2 are (i = 1,2) 

(Hi - ^jCz^ Ui - giiljfvi + gi2'ipin'^-i = hcjUi, (2) 
(Hi + QCz^ Vi - gi^pfui - gi2'^*n'r^-i 



-hcovi, 



where Hi = H - j^i ^ '^9i\'4^i? + ^i2|V^3-iP with i = 
1,2 and n[ = — il^iVi is the density fluctuation 

for species i p^. We use a harmonic trap potential 
Vtr = MujI^{x^ -h ?/2)/2 + Mulz'^/2. We further as- 
sume Qi = g2 = g which is relevant for the experi- 
ment 9]. The BdG equations are solved numerically by 
discretizing the eigenvalue problem, using finite element 
method with local cubic Hermite polynomials as the ba- 
sis. The solutions are analyzed by observing the dynam- 
ics of ilJi{t) = il)^ + UiC'^"^^ - v*e*'^^ 11] together with 
the density fluctuation n-. Here i/j^ is the ground state 
solution of the i-th component. Our analysis assumes a 
quasi 2-dimensional geometry appropriate for a pancake 
shaped trapping potential with the dynamics along the 
z-axis frozen. 

With gi2 < g^ the energy spectrum of a non-rotating 
PO condensate is separated into two branches with in- 
phase (+) and out-of-phase (— ) modes. Each mode 
can be characterized using radial and angular quantum 
numbers, n and m just as is done for a single com- 
ponent condensate O Q 10- The overall density 
rit = V^iV^i + '02'^2 for the out-of-phase modes does not 
oscillate significantly, which can be understood by noting 
that each species moves oppositely to the other, remi- 
niscent of Landau's second sound in superfluid helium. 
To concretely distinguish the (+) modes from the (— ) 
modes, the difference in the complex phase between two 
density fluctuations n[ and n2 can be quantified by defin- 
ing D = J dYnin'2/ J dr\n[\\n2\ which takes the value 
between +1 and —1:D ~ 1 for (+) modes and —1 for 
(— ) modes. Figure shows the departure of the (— ) 
modes (blue crosses) from (+) modes (red dots) as gi2 
increases from zero up to gi2 = g. The frequencies of the 
(+) modes do not change significantly with the change 
of gi2. To interpret this, note that the (+) modes are 
barely affected by the interspecies interaction as the two 
species move together rather than against each other. 
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FIG. 1: Red dots and blue crosses represent in-phase and out- 
of phase modes respectively. Out-of-phase modes gradually de- 
part from in-phase modes decreasing energies as interspecies 
interaction increases. The spectrum of in-phase modes are not 
seriously affected by change of gi2- 




FIG. 2: Three ground states depending on different rotational 
frequencies with gi2 — 0.75g. Bright clouds represent one 
species with red dots the positions of vortices while blue dots 
are vortices of the other component. 



As in a single species case [13 it is found that highly 
excited surface modes also exist that involve a circle of 
phase singularities flowing together with a density wave 
around the edge of the condensate. For the (+) modes 
the phase singularities for each of the two species over- 
lap one another, while for the (— ) modes, they alternate 
around the edge of the condensate. For a given quantum 
number m, the (— ) mode has lower energy compared to 
the (+) mode since, for gi2 > 0, the separation between 
phase singularities from different species is energetically 
preferred. 

When a PO BEG is rotated, it develops two interlacing 
rectangular vortex lattices - one condensate fills up the 
density holes made by the vortices in the other compo- 
nent [Fig. 12]. It is interesting to note that, because of the 
finite number of vortices corresponding to a given rota- 
tional frequency, the azimuthal symmetry of the vortex 
lattice structure differs for different rotational frequen- 
cies. As in a rotating single component BEG, two types 
of normal modes were found to emerge once the system 
contains vortices: the excitations of the vortex structure 
(Tkachenko modes) and fluid density fluctuations of the 
underlying condensate (hydrodynamic modes) including 
surface modes. The surface modes constitute one of the 
main types of excitations in the lower energy regime as 
shown in Fig. [31 the energy vs. angular momentum plots 
for gi2 = 0.75g. For the Tkachenko modes, one might ex- 
pect that, with two species present, there would be two 
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FIG. 3: Energy vs. angular momentum for the three 
ground states in Fig. |21 Each point represents a normal 
mode. Tkachenko modes are identified as negative angular mo- 
mentum states in the low energy regime. Surface modes 
are divided into in-phase (red dots) and out-of-phase (blue 
crosses) modes. Note that the gap between spectral lines of 
in-phase and out-of-phase surface modes decreases at higher 
angular frequencies. 



branches of normal modes, the "optical" and "acoustic" 
modes similarly to a solid state crystal lattice with two 
constituents [la^llj]- Numerical calculation, however, in- 
dicates that the distinction between these two branches is 
not evident. This is due to the fact that in the BEG sys- 
tem the interaction between vortices in different species 
is caused by density- density interactions and hence is far 
smaller than the inter- vortex interaction within the same 
species dominated by the kinetic energy term. 

Figure El shows that the discrepancy between in-phase 
(red dots) and out-of-phase (blue crosses) modes becomes 
negligible as ft increases. This is because the phase sin- 
gularities involved in the surface modes no longer have 
either clean overlap or separation due to the presence of 
vortices near the edge of the condensate. When the vor- 
tex lattice structure near the edge of the condensate is 
azimuthally asymmetric, the surface wave of one species 
interweaves the surface wave of the other species. For 
the case of azimuthally symmetric vortex lattice struc- 
ture, however, the phase singularities that constitute the 
surface wave of one species is uniformly closer to the edge 
of the condensate than those of the other species due to 
the repulsive interaction between the outer vortices of the 
vortex lattice and the phase singularities (vortices) in the 
surface wave. It was found that in such cases, the ratio of 
the quasiparticle populations from the two species inte- 
grated through the whole condensate departs from unity 
dramatically resulting in only one of the species fluctuat- 
ing while the other component remains nearly stationary: 
the quasiparticle population of the species whose excita- 
tion is closer to the edge of the condensate dominates 
the other. Such quasiparticle population imbalance of 
the surface modes are a new feature of the rotating con- 
densate, and can provide a potential diagnostic tool for 
the vortex lattice structure. 

For gi2 > g, due to the broken rotational symmetry 
of the non-rotating PS BEG ground state, the degenera- 
cies of +m and — m states of the PO {gi2 < g) state for 
the out-of-phase modes no longer persists. It has been 
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FIG. 4: Energies of normal modes as a function of gi2 show- 
ing the transition from phase-overlapped to phase-separated 
state. The actual transition interaction strength is slightly 
larger than g due to the kinetic energy contribution. 



suggested that due to the surface tension there might 
therefore be a new branch of normal modes on the sur- 
face that separates the two species 6]. Our numerical 
solution shows that, for a strongly phase-separated con- 
densate, such modes exist and indeed dominate the lower 
energy regime. These modes are not supplemental to the 
hydrodynamic excitations; instead they can be extrapo- 
lated from the out-of-phase hydrodynamic modes of PO 
BEG by changing gi2 continuously. Figure ^ shows the 
change in spectrum as the interspecies interaction goes 
through the transition point. The actual transition point 
is slightly larger than g due to the kinetic energy con- 
tribution. Note that, for gi2 > g, there are different 
slopes for different branches, as in a single-species BEG 
with broken rotational symmetry . The group of lines 
with the lowest slope (~ 0) corresponds to the intersur- 
face wave propagating between species. They have differ- 
ent numbers of nodes depending on the intercept of the 
spectral lines with energy axis. The frequencies of these 
modes do not noticeably depend on any form of surface 
tension (or ^12) as they do not penetrate far into the 
condensate. The modes with non-zero slopes i.e. those 
with energies which depend on gi2 have additional nodes 
in the density fluctuation along the axis perpendicular 
to the inter-surface and the slope is determined by the 
number of such nodes. 

When the PS condensate is rotated, the formation of 
vortex sheets in the system effectively slices the conden- 
sate into several layers. By forming vortex sheets, the 
condensate achieves maximal phase separation compared 
to a lattice. It is notable that, for a given vortex den- 
sity (which depends on the mass of the atom and the 
rotational speed of the system) one can achieve much 
smaller intervortex spacing within a vortex sheet when 
compared to a lattice geometry. We find that, for the 
case of vortex sheets, in contrast to the unique structure 
of a vortex lattice, there exist different met ast able vortex 
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(c) Q = O.ScJt: 



FIG. 5: Three ground states depending on different rota- 
tional frequencies with gi2 = l-ls'- Bright clouds represent 
one species with red dots the positions of vortices while blue 
dots are vortices of the other component. (a) is close to the 
nonrotating PS condensate except that it has phase singular- 
ities stabilized by the presence of the other component. (b) is 
the first emergence of a vortex sheet (red dots) that effectively 
divides the condensate into two pieces, (c) is a typical vortex 
sheet state for high rotational frequency. 
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FIG. 6: Spectrum of three different vortex sheet states corre- 
sponding to Figs. |5|(a-c). 



sheet configurations with almost the same energies. We 
understand this by noting that the energy of the system 
is mainly determined by two length scales, the intervor- 
tex spacing (dy) within a vortex sheet and the intersheet 
spacing (d^), rather than the shape of the vortex sheets 
However, since any bending or defects in the sheet 
configuration still has some energy cost, the state with 
straight vortex sheets corresponds to the lowest energy 
for our parameters and has been used for the analysis 
[Fig.ISl. 

Figure shows the frequency vs. angular momen- 
tum curve for gi2 = l.lg. The Tkachenko modes are 
clearly identified in the lower energy regime especially 
for the fast rotating case (c), while hydrodynamic ex- 
citations are found throughout the whole region of en- 
ergy. It is obvious that the wavelength of the Tkachenko 
mode should be bounded below by dy. We find that, 
different from the triangular or rectangular vortex lat- 
tices, Tkachenko modes are anisotropic due to the in- 
homogeneity of the vortex matter. The resulting den- 
sity fluctuations are now concentrated along the vortex 
sheet. Figure [3 shows |n'p for two Tkachenko modes, 
with wavelengths longer and shorter than the dg respec- 
tively. The long-wavelength Tkachenko mode [Fig. HJa)] 
is excited homogeneously throughout the whole region of 
the vortex matter since the wavelength is longer than 
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FIG. 7: Amplitude of density fluctuation n of Tkachenko 
modes with wave length (a) longer and (b) shorter than in- 
tersheet spacing. 




FIG. 8: A snapshot of the most highly excited Tkachenko 
mode. Arrows are added to show the motion of vortices flow- 
ing along each sheet. 

(is, while the short- wavelength Tkachenko mode [Fig. 
El^b)] is excited mainly along a single sheet with the most 
anisotropy since waves with wavelength shorter than ds 
cannot travel across from one sheet to another. 

We have also observed the breakdown of stable vor- 
tex positions associated with the excitations of certain 
Tkachenko modes. This happens when the typical vortex 
precession radius becomes the same order as dy which, 



in turn, is comparable to the healing length of a vortex 
within a vortex sheet. For some highly excited Tkachenko 
modes, the relative precession motion of the nearest 
neighbor vortices are out of phase such that the vortices 
"collide" with one another when the precession radius is 
larger than dy/2. The result is an exchange of vortex 
positions in the lattice site. The successive exchange of 
vortex positions leads to a shear flow of "delocalized" 
vortices along the vortex sheet [Fig. [Sj. We have con- 
firmed that this flow can happen stably and without ex- 
periencing any significant nonlinear decay within a time- 
dependent GPE simulation. This vortex flow gives rise to 
"inner" surface waves along the edge of the layers of the 
condensate and therefore such modes may provide a con- 
nection between surface modes and Tkachenko modes. In 
general these "inner" surface waves cannot flow out nat- 
urally to the outer surface of the overall BEG since the 
frequencies and wavelengths of the "inner" and "outer" 
surface waves do not necessarily match. However, for 
some modes, we have observed that the two types of sur- 
face waves can be connected through the generation of 
temporary vortex- ant ivortex pair which provides the nec- 
essary frequency matching. 

In summary, we have studied the normal mode dynam- 
ics of a rotating two-species BEG and found new features 
that result from the different ground state structures. For 
PC state with rectangular lattices, we found that quasi- 
particle population imbalance may provides a diagnostic 
tool for the vortex lattice structure. For PS state with 
vortex sheets, we observed the breakdown of ground state 
vortex positions through shear flow of vortices within the 
linearized regime. We note that as compared to a vortex 
lattice with similar vortex density, the shorter intervor- 
tex distance in the sheets allows observation of vortex 
collisions at lower rotation frequency. The potential link 
between vortex delocalization and vortex lattice melting 
will be a topic of future investigation. 
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